Nonresonant B\ g Raman scattering in the Hubbard model 



O 

o 

(N 
u 

< 



■4— » 

c3 



c 

o 
o 



> ■ 

00 ' 

^f- : 
o . 

o ■ 



E 

T3 ■ 

c : 
o . 
o ■ 



X 



J.K. Freericks*, T.P. DevereauxT and R. Bulla* 

* Department of Physics, Georgetown University, Washington, DC 20057, U.S.A. 
^Department of Physics, University of Waterloo, Canada 
* Theoretische Physik HI, Elektronische Korrelationen und Magnetismus, Institut fur Physik, Universitdt Augsburg, D-86135 

Augsburg, Germany 
(February 1, 2008) 

The numerically exact solution for nonresonant B\ s Raman scattering is presented for the half- 
filled Hubbard model in infinite dimensions. This solution illustrates the modifications of the Raman 
response (in a system tuned through the quantum-critical point of a metal-insulator transition) due 
to Fermi-liquid properties in the metallic phase. In the insulating phase, we recover the predicted 
universal behavior, while we find the Raman response is quite anomalous on the metallic side of 
the transition. Our calculated results are similar to those measured in FeSi, SmB6, and underdoped 
cuprates. 



Principle: 78.30.-j; 71.30.+h; 74.72.-h; 75.20HR 

Raman scattering has been applied to study the excita- 
tions of electrons in metals, insulators, semiconductors, 
and superconductors. Inelastic light scattering reveals 
symmetry selective properties of the electron dynamics 
over a wide range of energy scales and temperatures (the 
symmetry is selected by polarizing the incident and re- 
flected light). Three principle channels are usually exam- 
ined: (i) A\ g which has the full symmetry of the lattice 
(i.e. is s-like); (ii) B\ g which is a d-like symmetry (that 
probes the Brillouin-zone axes); and (iii) £?2g which is 
another d-like symmetry (that probes the Brillouin-zone 
diagonals). Shastry and Shraiman (SS) proposed a sim- 
ple relation [Q that connects the nonresonant Raman re- 
sponse to the optical conductivity, which was recently 
proven in infinite-dimensions Q for the B\ g channel. We 
employ the SS relation here to evaluate Raman scattering 
in the Hubbard model. 

Strongly correlated systems as disparate as mixed- 
valence compounds || (such as SmB 6 ), Kondo- 
insulators [Q (such as FeSi), and the underdoped 
cuprate high temperature superconductors show 
temperature-dependent Bi g Raman spectra that are 
both remarkably similar and quite anomalous (see Fig- 
ure [l]) . This "universality" suggests that there is a com- 
mon mechanism governing the electronic transport in 
correlated insulators. As these materials are cooled, they 
all show an increase in the spectral weight in a "charge- 
transfer" peak with a simultaneous reduction of low- 
frequency spectral weight. This spectral weight trans- 
fer is slow at high temperatures and then rapidly in- 
creases as the temperature is lowered in the proximity 
of a quantum-critical point (corresponding to a metal- 
insulator transition). The Raman spectral range is also 
separated into two regions: one where the response de- 
creases as T is lowered and one where the response in- 
creases. The characteristic frequency that divides these 
two regions is called the isosbestic point, which is the 
frequency where the Raman response is independent of 
temperature. These anomalous features are not typically 



seen in either the Ai g or the i?2g channels. 

This anomalous experimental behavior is shown in Fig- 
ure [j]. The top panel shows SmBg, which has the added 
feature of developing a sharp peak at 130 cm -1 (that 
does not disperse in frequency) when the temperature is 
lower than 30 K. The FeSi data is shown in the middle 
panel. It displays the cleanest signature of these anoma- 
lous features. Note how the isosbestic point only develops 
at temperatures below 150 K. The bottom panel shows 
the data in the LSCO high-temperature superconductor 
(which we have smoothed using a 4253H, twice smoother). 
The isosbestic point is somewhat harder to see here (be- 
cause of the noise in the data) , but it does clearly develop 
at about 2000 cm -1 as the temperature is lowered (one 
cannot go below the superconducting transition temper- 
ature, since additional signatures of the ordered phase 
would be seen in the Raman response). 

Theory has just begun to catch up with experiment for 
electronic Raman scattering in strongly correlated mate- 
rials that are tuned to lie near a metal-insulator transi- 
tion This theoretical treatment involved applying 
the dynamical mean-field theory to the simplest many- 
body system that has a metal-insulator transition — the 
spinless Falicov-Kimball model || . While results for the 
Big channel were shown to be universal on the insulat- 
ing side of the transition (and qualitatively agree with 
the anomalies seen in experiment), on the metallic side, 
the effects of a low-temperature Fermi coherence peak are 
unknown, since the Falicov-Kimball model is a non Fermi 
liquid. The iterated-perturbation-theory approximation 
has been applied to the Hubbard model before |[o| but 
not in the regime of the metal-insulator transition; here 
we provide the numerically exact solution in this regime. 

The nonresonant Raman scattering in the £?2 g channel 
vanishes and in the Ai g channel it requires knowledge 
of the local irreducible charge vertex (which is problem- 
atic to calculate in the Hubbard model due to its large 
dynamic range, which cannot be treated by maximum- 
entropy techniques), so we provide results only for the 
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FIG. 1. Experimental Raman response for correlated 
materials (a) SmBg [3]; (b) FeSi [4]; and (c) underdoped 
La2- a; Sr :E Cu04 [5] with x = 0.08. All of the experimental 
data show the development of a low-temperature isosbestic 
point, which occurs due to the transfer of spectral weight 
from low energy to high energy as the temperature is low- 
ered, indicating the proximity to the quantum-critical point 
of a metal-insulator transition. The individual curves are la- 
beled by the temperature in K where the measurement was 
taken. In panel (c) only the high temperature (300 K) and 
the low temperature (50 K) are labeled. The two intermediate 
curves are at 100 and 200 K respectively. 



B\ g sector here. The Hubbard Hamiltonian contains two 
terms: the electrons can hop between nearest neighbors 
[with hopping integral t* /(2yd) on a d-dimensional hy- 
percubic lattice p!l|], and they interact via a screened 
Coulomb interaction U when they sit on the same site. 
All energies are measured in units of t* . The Hamiltonian 
is 



H = - 



— T 



C ia C 3° 



(1) 



where cj CT (ci a ) is the creation (annihilation) operator for 
an electron at lattice site i with spin a and = cl a Ci a 
is the electron number operator. We adjust a chemical 
potential fi to fix the average filling of the electrons. 

The nonresonant Raman response in the Bi g channel 
has no vertex corrections []l2|,[s|]2"[ and is equal to the bare 
bubble. The formula for the Raman response has been 
presented elsewhere Q and can also be found from the 
SS relation Jj],[l3j — the imaginary part of the nonresonant 
Bi g Raman response is 

1mR(u) =c J duj[f(uj) - f(u> + v)] 

x / dep(e)A(e,oj)A(e,uj + v), (2) 



where f(u>) = 1/[1 + exp(w/T)] is the Fermi function, 
A(e,oj) — Im[— l/ir{ui + /i — £(o>) — e}] is the spectral 
function, p(e) is the noninteracting density of states (a 
Gaussian here), and c is a constant. 

We study the evolution of the Raman response at 
half filling, since one can tune the system to move right 
through the quantum-critical point of the metal-insulator 
transition. We analyze finite-temperature numerical 
renormalization group calculations []l4| (restricted to the 
paramagnetic phase which has a metal-insulator tran- 
sition). We examine three cases here: (i) a correlated 
insulator just above the transition U — 4.24 (where the 
response is "universal"); (ii) a metal just below the phase 
transition U = 3.54 (which undergoes a temperature- 
dependent metal-insulator transition at 0.011); and 
(iii) a correlated metal U — 2.12. 

We show the correlated insulator regime in Fig. ||, 
where the Raman response is model-independent. We 
see behavior identical to that seen in the Falicov-Kimball 
model (the only modification here is the consequence 
of a temperature dependence in the interacting density 
of states, which fills in the gap at temperatures above 
about 0.1). We see the development of low-energy spec- 
tral weight at the expense of the higher-energy charge- 
transfer peak as T increases (although in this case the 
spectral weight grows over a broader temperature scale) 
and the appearance of a single isosbestic point at a some- 
what higher energy than seen in the Falicov-Kimball 
model solution [y « U/1.5 rather than v « U/2). All 
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FIG. 2. Nonresonant Bi g Raman scattering for a correlated 
insulator U = 4.24 and half filling for a number of tempera- 
tures ranging from 0.57 to 0.071. Inset is an enlargement of 
the low-energy features. 
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FIG. 3. Nonresonant Bi s Raman scattering for system 
undergoing a temperature-driven metal-insulator transition 
U = 3.54 and half filling for a number of temperatures ranging 
from 0.28 to 0.003. Inset is an enlargement of the low-energy 
features. 



of the qualitative features of the universal behavior are 
shared in the Hubbard-model solution. The difference 
is in the quantitative details, the most apparent one be- 
ing that at high temperatures the low-energy spectral re- 
sponse does not have a broad peak structure, but rather 
shows monotonic rising behavior. This is more charac- 
teristic of experiments in materials like FeSi Q or the 
cuprates which have a relatively flat low-energy 

spectral response. 

In Fig. ||, we show results for a system tuned to lie 
just on the metallic side of the metal-insulator transi- 
tion (so it undergoes a temperature-driven transition at 
T 0.011). Initially the Raman response acts like an in- 
sulator: the charge-transfer peak sharpens and grows in 
strength, while the low-energy weight is reduced and an 
isosbestic point appears near v = U/1.7 as T is lowered. 
As the temperature is reduced further, the response be- 
comes quite anomalous. Spectral weight shifts back out 
of the charge-transfer peak into the low-energy region, 
developing two low-energy bumps and an isosbestic point 
near v = 2. There is only one experimentally measured 
system that we know of that shares some of these qual- 
itative features: SmBg ||. At low temperatures SmBg 
acts like an insulator, just as seen in Figure ||, but then 
at the lowest temperatures (T < 30 K), a sharp nondis- 
persive peak appears at about 130 cm -1 . The qualita- 
tive shape of the Raman response is different though, 
because the weight only grows in a narrow peak, as op- 
posed to the wide frequency range of Figure ||. We be- 
lieve it would be interesting to investigate the Raman 
response of a correlated system that undergoes a simi- 
lar insulator-metal transition as a function of tempera- 
ture such as 1% Chromium doped Vanadium Oxide [ jT5| . 
Such experiments would be feasible if the charge-transfer 
peak could be pushed to a low enough energy that it lies 



within the window observable by Raman measurements 
(which seems likely from the optical conductivity data 
on insulating Vanadium Oxide alloys doped close to the 
metal- insulator transition). 

Finally, in Fig. ||, we show the Raman response for 
a correlated metal at half filling. At high temperatures, 
the Raman response has a wide charge-transfer peak cen- 
tered at u w U. As the temperature is lowered, we see 
the development and evolution of a low-energy Fermi- 
liquid peak, which sharpens as T is lowered. This is the 
classic behavior expected for a correlated metal — at high 
temperatures there is a large charge-transfer peak cen- 
tered at an energy just somewhat higher than U which 
loses spectral weight as the temperature is lowered and a 
low-temperature "metallic" peak develops at low energy. 
The Fermi peak has the expected form proportional to 
vaT/(v 2 + a 2 T 2 ) for a Fermi liquid. The width of the 
peak (determined by aT) decreases as the temperature is 
lowered and will ultimately vanish at T = 0. The weight 
in the metallic peak is much smaller than the weight in 
the charge-transfer peak. The continuous evolution of 
the low-temperature peak down to zero frequency is not 
seen in the Falicov-Kimball model, because the self en- 
ergy remains finite there in the limit as T — > 0. Note 
that the charge-transfer peak loses weight (especially on 
the low-frequency side) as T is lowered, but there is no 
low-energy isosestic point here (in fact an isosbestic point 
may be developing at v w 3.3). This lack of low-energy 
isosbestic behavior is quite interesting. The development 
of a large- weight Fermi coherence peak in the interacting 
density of states tends to destroy the isosbestic behavior. 

Unfortunately, we are not aware of any experimen- 
tal data on a correlated metal that displays this low- 
temperature development and evolution of the Fermi- 
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FIG. 4. Nonresonant Bi g Raman scattering for a correlated 
metal U = 2.12 at half filling for a number of temperatures 
ranging from 0.35 to 0.024. Inset is an enlargement of the 
low-energy features. 

liquid peak. Surprisingly, little is known about Raman 
scattering in a weakly interacting metal. It is well know 
that in the absence of inelastic scattering, the low energy 
Raman cross-section vanishes at q=0 due to the lack of 
phase space available to create particle- hole pairs. There- 
fore any signal at all must come from electron-electron 
interactions fl6|| or impurities |l7| . Our results show for 
the first time that the consequences of well-defined quasi- 
particles is the presence of a low energy peak which grows 
in intensity and narrows as temperature decreases. Such 
a peak might also be observable in materials like CeSi2, 
CeSb 2 , CeBei 3 (l| or YbAl 3 (l| which all display the 
development of a low-temperature Fermi coherence peak 
at energies less than 100 meV, but may be below the 
sensitivity of most experiments near the laser line. 

To summarize, we have examined the effects of Fermi- 
liquid behavior in the metallic phase on the Raman re- 
sponse of a correlated material that is tuned to lie close 
to the quantum-critical point of a metal-insulator tran- 
sition. In the case of half-filling, where we can tune the 
system to move right through the quantum-critical point, 
we see similar behavior as seen in the Falicov-Kimball 
model, but influenced by the appearance of Fermi-liquid 
behavior on the metallic side of the transition. Hence, for 
a correlated metal, there is a charge-transfer peak, which 
has a Fermi-liquid peak separate from it at low tempera- 
tures and frequencies. The peak sharpens and moves to 
lower energy as T is lowered and an isosbestic point may 
be forming on the high-energy-side of the charge-transfer 
peak. In the insulating regime, we recover the universal 
behavior predicted, which shows a low-temperature in- 
crease in low frequency spectral weight at the expense of 
the charge-transfer peak and the appearance of a single 
isosbestic point (but with somewhat stronger tempera- 
ture dependence) . The quantitative features appear to be 



closer to experiment than what was seen in the Falicov- 
Kimball model. The region closest to the quantum- 
critical point is the most interesting. In the insulating 
(high-temperature) phase it behaves like a correlated in- 
sulator, but once it develops a low temperature Fermi- 
liquid peak in the density of states, the low frequency 
Raman response increases as T is lowered, and develops 
a series of low-energy bumps. This behavior hasn't yet 
been seen experimentally, but it would be interesting to 
look for in a material that undergoes a low temperature 
metal-insulator transition (such as the doped Vanadium- 
Oxide system). 
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